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Abstract: We study the continuum limit of the string- hke behaviour of flux tubes formed 
between static quarks and anti-quarks in three dimensional SU{2) lattice gauge theory. We 
compare our simulation data with the predictions of both effective string models as well 
as perturbation theory. On the string side we obtain clear evidence for convergence of 
data to predictions of Nambu-Goto theory. We comment on the scales at which the static 
potential starts departing from one loop perturbation theory and then again being well 
described by effective string theories. We also estimate the leading corrections to the one- 
loop perturbative potential as well as the Nambu-Goto effective string. In the intermediate 
regions we find that a modified Lennard-Jones type potential gives surprisingly good fits. 
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1. Introduction 

Gluonic dynamics at large distances can be very effectively probed by looking at proper- 
ties of the flux tube which forms between a static quark and an anti-quark in the QCD 
vacuum. Although there is still no analytical proof of this formation, lattice simulation 
results overwhelmingly indicate that this indeed is the case Q. 

Bosonic string descriptions of this flux tube have been around for a long time. While 
most of the earlier attempts Q used an open string description with the ends of the string 
ending on the quark and the anti-quark, recently there have been attempts to give closed 
string descriptions too. After Polchinski and Strominger(PS) Q suggested how effective 
theory of strings with vanishing conformal anomaly could be formulated in all dimensions 
it has been shown |5|, ^ that the spectrum of these effective theories is universal^ to 
order (r being the length of the string) and that to this order they coincide with the 
predictions of Nambu-Goto theory whose conformal anomaly however vanishes only in 26 
dimensions. 

Another interesting idea put forward by Liischer and Weisz 0] is that of open-closed 
string duality. They showed that this too constrained the possible string spectra. In three 
dimensions they found that this duality implied that to order the spectrum was the 
same as that of the Nambu-Goto(NG) string while in four dimensions one parameter was 
left undetermined. In the PS effective string theories the spectrum to this order is the same 
as NG theory in all dimensions without needing to invoke open-closed duality explicitly. 

On the lattice, the flux tube can be observed as the potential between a static quark 
and an anti-quark (open string). This can be obtained from expectation values of Wilson 
loops or Polyakov loop correlators. One of the characteristics of the string like behaviour 
of the flux tube is the presence of a long distance 1/r term in all dimensions and with 
a universal coefficient in the qq potential. In four dimensions there is a short distance 

^In this article universal will mean independent of the details of the underlying gauge group. 
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l/rlog(r) term in the potential. However its coefficient depends on the details of the gauge 
group distinguishing it from the long distance 1/r term. The long distance 1/r term was 
first observed in ^ and its universality was established in [^]. It is known as the Liischer 
term. 



The Liischer term has been looked at in lattice simulations since the eighties |1C, 11|. 
However in recent times, increase in computing power and improvement in algorithms have 
allowed really precise measurements of that term and also the term. Now it is possible 
to study the properties of flux tubes longer than 1 fermi which was unimaginable even a few 
years back and one is really in a position to compare the lattice data with the predictions 



of the string picture. See ||12|, |14|, |16|, [l^, |19|, for example. In particular, 



measurements of the Liischer term in three dimensions have been carried out in |0, 19 1 for 
SU{2) , in |14| for SU{3) and in [|^] for SU{5) lattice gauge theories. Convergence to the 
static potential of NG theory to order at a distance scale of around 1 fermi has been 



recently shown in the case of d = 4 SU{3) gauge theories in |15|. 

Another characteristic of the string behaviour of the flux tube is the level spacing of 
the excitation spectrum. We do not discuss it here. See 0, |l^, 17, 21, 22| for recent analytic 



and numerical studies. A review on this topic can also be found in |23|. 

In this article we present results of our simulations of the Polyakov loop correlators 
for d = 3 SU{2) Yang-Mills theory and compare the resulting static potential with both 
perturbation theory and string model predictions. This allows us to narrow down bounds 
on the distance beyond which we can say the flux tube indeed shows a string like behaviour. 

2. Simulation parameters 

We have carried out simulations of three dimensional SU (2) lattice gauge theory on lattices 
at four different lattice spacings with the coarsest lattice having a spacing of slightly below 
0.13 fm to the finest lattice with a spacing of about 0.045 fermi. We use symmetric cubic 
lattices and the Wilson gauge action. On all these lattices, we have computed Polyakov 
loop correlators {P*{x)P{y)) for various spatial separations r = y — x. 

To reliably extract signals of these observables which are exponentially decreasing 
functions of r and T (the temporal extent of the lattice), we used the Liischer- Weisz 



exponential error reduction algorithm |13]. In this algorithm, one computes intermediate 
expectation values on sub-lattices of the original lattice which are obtained by imposing 
suitable boundary conditions. For measuring Polyakov loop correlators, we obtain our 
sub-lattices by slicing the original lattice along the temporal direction. As is well known, 
this algorithm has several optimization parameters among which the number of sub-lattice 
updates employed seems to be the most important one. The lattice parameters along with 
the number of sub-lattice updates used in each set of measurements are summarized in 
table |. 

Another important parameter is the thickness of the time-slice over which the sub- 
lattice averages are carried out. We found that it was helpful to increase this thickness as 
one goes from stronger to weaker coupling. We used values of two, four and six as time-slice 
thicknesses. 



- 2 - 



(3 r values lattice iupd # of measurements 



5.0 


2-8 


36^ 


16000 


1600 




7-9 


40^ 


32000 


3200 




8-12 


48'^ 


48000 


7000 


7.5 


4-8 


48'^ 
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7-12 


643 


18000 


1100 




11 - 16 


64^ 


36000 


7200 


lU.U 


2 — 7 


48 


16000 


2850 




6-9 


48^ 


16000 


200 




8-14 


843 


24000 


1100 




13-19 


843 


36000 


2250 


12.5 


2-9 


48^ 


16000 


2700 




8-14 


723 


24000 


1150 



Table 1: Runs 



To optimize the running time while keeping any finite volume effect under control we 
chose different lattice sizes for different ranges of r. Typically the smaller values of r are 
much easier to get as they can be reliably obtained on smaller lattices and require much 
less sub-lattice updates. In principle of course results for all the different r values could be 
obtained from the largest lattice but memory requirements prevent us from doing all the 
measurements in a single run. 
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3.9536 (3) 
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1.5214 (3) 
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0.861665 (4) 
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6.2875 (10) 


0.07952 


(1) 


1.5246 (7) 


10.0 


0.897683 (3) 


0.020606 (4) 


8.6022 (8) 


0.05812 


(1) 


1.5248 (4) 


12.5 


0.921100 (2) 


0.012742 (17) 


10.916 (3) 


0.04580 


(1) 


1.5215 (29) 






Table 2: 


Lattice scales 









The scale is set by the Sommer scale tq = 0.5 fm which is implicitly defined by 
rQ/(ro) = 1.65, where /(r) is the force between the static quark and the anti-quark. We 
use our measured force values and interpolation to extract tq. The values that we obtain 
are given in table ^. We also see from table || that otq is a constant (~ 1.522) to a very 
good approximation as expected in the scaling region close to the continuum limit. 

3. Results 

In this section we present the results of our measurements. From the {P*P) correlator one 
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can extract the static quark-antiquark potential V{r) by 

y(r) = -iln(P*P(r)). (3.1) 

In principle the qq potential contains all the information about the flux tube, but it also 
contains an unphysical constant. We therefore look directly at the first and the second 
derivative of this potential. 

The first derivative of the potential gives us the force between the quark and the anti- 
quark, while the second derivative gives us information about the subleading terms and how 
one approaches the asymptotic linearly rising behaviour of the qq potential. To facilitate 
our comparison with string models we actually compute a scaled second derivative which 
we call c(f). This quantity should become the Liischer term (= —-^^21^) asymptotically. 

On the lattice these quantities are given by 

F{r) = V{r) - V{r - 1) (3.2) 

~3 

<f) = y [y (r + 1) + V{r - 1) - 2V{r)] (3.3) 

where f = r + | + O(a^) and f = r + O(a^) are defined as in [0] to reduce lattice artifacts. 

The theoretical predictions in continuum are given by string models (for large r) as well 
as perturbation theory (for small r). Since non-bosonic string models have been essentially 



ruled out |2j], we are going to concentrate on the potential due to the NG string, the so 



called Arvis potential |25] given by 



(d-2)7rV/2 



y^r.. = ar^l-^^) . (3.4) 

Keeping in mind the results of Liischer and Weisz ||^ and PS type effective string theories 
1^, ^, ^ we are going to compare our lattice data on force and c(r) with leading order 
predictions to which all models reduce at sufficiently large r, to NLO expressions and 
expressions from the full Arvis potential. In three dimensions the expressions are given by 



L.O. /(r) = <r+(^)i c(r) 

A™s /(.) = . (l - c(r) = (3.5) 
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The perturbative potential has been calculated at the one loop level by Schroder pq ]. 
He obtains 

Vpertif) = Cpcrt?" H In g^r + (higher order terms) (3.6) 

with the perturbative string tension cXpert = '^^~^r^- For SU{2) Cp = 3/4, Ca = 2. 
The perturbative force and Cpert(?') can be computed by /pert(^) = '^^''^^ ^"^^ Cpcrt('') = 

r3 d?Vpcrt{T) 

2 dp- • 



- 4 - 




1 I < < < < 1 

0.5 1 1.5 2 2.5 3 

r/ro 

Figure 1: r^fir) vs r/ro in the d = 3 SU{2) case. The four different sets are D {(3 — 5), 
V (/3 = 7.5), A (/3 = 10) and O {(3 = 12.5). Also shown is the 1-loop perturbation theory curve 
/port(r), as weh as the leading order (L.O. /(r)) and NG (Arvis /(r)) string predictions. The 
horizontal line is ro/(r) ~ 1.65 and locates the Sommer scale. 



We first look at the force data. Since we have four values of the coupling, we look at 
the continuum limit of the string tension. Following ||2^ we too define /3mf = /? x ^(Tr up) 
and look at how a^/a scales with /?mf- We use /3mf = ^/g^a to convert the perturbative 
expressions to lattice units. Fitting the string tension to the form 



where N refers to the gauge group SU{N), we obtain (^) , for iV = 2, to be 

V^' / cont. 

0.16788(12) which is completely consistent with the values presented in [27, 25]. As noted 
by Schroder the full string tension is 1.47 times the perturbative string tension and 
hardly changes for between 2 and 5. We want to also mention that if we ignore the 
point /3 = 5, then the data can be well described even without the coefficient a2- In that 
case we obtain the continuum value to be 0.16736(10). Another point worth noting is that 



while we too find the coefficient ai to be negative, in agreement with |27], the coefficient 
02 is positive in our case. This is consistent with the higher A^'s in [^^, and may not be of 
much significance as in that work, 02 for SU (2) is consistent with within 2a. 

The values obtained for otq by fitting the force to the Arvis form is quoted in table |2|. 
The values given by the other two forms differs by less than 0.1% from the quoted values. 
ar^ scales very nicely with the different beta values differing by less than 0.2% from each 
other. The fits were carried out in the range 2.4 - 2.9 rg for /3 = 5 , 2.1 - 2.5 rg for (3 = 7.5, 
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1.9 - 2.1 ro for /? = 10 and finally between 1.6 - 1.8 tq for [i = 12.5. The data for large r 
at P = 12.5 was taken from fl% . 

In Fig. H we plot Tq/ versus r/rg which is expected to become a universal curve in 
the continuum limit. The horizontal line is r^/(r) = 1.65 and defines the Sommer scale 
tq. The grey line is the L.O. result otq + ^ (^) • The black line is the Arvis curve 
and the magenta one the perturbative curve. The data starts departing from the one-loop 
perturbative curve around 0.22 ro or 1.8 GeV and joins onto the string curves around 1.5 
ro or 260 MeV. The scaling exhibited by the data is very good with all the four different 
beta values falling on the same curve. This is mainly due to use of r instead of r , as it 
eliminates lattice artefacts to a large extent. The value for ar^ for drawing the L.O. and 
Arvis curves is taken from the fit at /? = 5 as it extends to the largest distance. Beyond 
1.5ro it is virtually impossible to distinguish the different theoretical curves as they are 
all dominated by the universal leading order behaviour with the string tension making up 
more than 95% of the force. The force data in fact gives the wrong impression of the string 
description being good even at distances as small as ro- 




Figure 2: c(f) in d = 3 SU{2) case. The four different sets are □ (/3 = 5), V (/? = 7.5), A {(3 = 10) 
and O — 12.5). Also shown are 1-loop perturbation theory {cpert{r)), L.O., N.L.O. and NG 
(Arvis) string model predictions. 

On the other hand, c(f) does not contain the string tension, and has a universal value 
in the L.O. It is therefore more sensitive to the sub-leading behaviour of the flux tube. In 
Fig. ^ we plot c(f) (versus r in units of ro) for all four /? values along with the perturbative 
curves as well as the three string model predictions given in eqn. The line and the 

symbol □ in red corresponds to /3 = 5, while the symbols V) ^ O colours green, blue 
and magenta correspond to the /? values 7.5, 10 and 12.5 respectively. The blue horizontal 
line corresponds to the leading order prediction while the ash and black curves correspond 
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to the NLO and Arvis forms respectively. We look at a wide range of r starting from where 
the data almost touches the perturbative curves going all the way to the region where the 
string predictions hold. 

The data seems to lie on top of each other exhibiting nice scaling behaviour as one 
goes to larger values of r. The /3 = 12.5 and /? = 10 data lie on top of each other in the 
range 0.75 and 1.25 tq already. The f3 = 7.5 set joins onto this at around 1.5ro and even 
the /3 = 5 data joins up at around 2.25ro. This points to the possibility that the continuum 
limit of the scale where the flux tube is well described by the Arvis curve can be obtained 
already on relatively coarse lattices. 



4. Discussions 



In d = 3, as also noted by Schroder, infrared divergences prevent computation of the 
perturbative potential beyond one loop. Infrared counterterms have to be obtained non- 
perturbatively As our data goes down to distances of about 0.16 fermi, we can try to 
obtain these counter terms by looking at two loop terms in the perturbative potential. On 
dimensional grounds, one can expect the perturbative potential to be of the form 

^pert(r) = ^ In g'r + ^^f^r + Ag'r Ing'r + Bg'r' + • • • (4.1) 
Ztt bill 

This gives 

cp=„{.) = + ^ (^)% i-t (4.2) 

47r ro 2 \rQj \roJ 

Since the first term is known, we determine A and B from the initial two points of our data 
to be A = 0.013162(3) and B = 0.001089(1). On our finest lattice we obtain g^ro, which is 
RG-invariant in the continuum limit, to be 3.792. While B contains effects of higher order 
terms, we expect the coefficient A to be relatively well determined. From the ratio of the 
leading to the next to leading term in Vpenir) as well as the data on c(f), we estimate the 
range of validity of first order perturbation theory to be about 0.1 fermi (consistent with 
our estimate from the force data). Fig. ^ suggests that second order perturbation theory 
holds upto distances of about 0.2 fermi. 

To try to get an idea about the scale of string formation, we look at the percentage of 
the total force carried by the string tension and the relative difference between the Arvis 
and the leading order force. From our data we find that the string tension constitutes 95% 
of the force at around 1.3ro, 98% at around 2.1ro and 99% at around 2.9ro. The relative 
difference, which gives us an idea about the importance of the subleading behaviour, is 
about 2% at 1.02ro, 1% at 1.2ro and goes down to 0.1% at about 1.9ro. 

The type of string is even more difficult to identify. At leading order, a variety of 



theories with different boundary conditions yield the universal Liischer term |29]. In fact 
all effective string theories of the PS type and AdS/CFT correspondences also yield 
this term. The type is therefore determined by the sub-leading behaviour of the flux tube. 
What can be clearly seen in the data is that the approach to the Liischer term is from 



^NDH wishes to thank G. 't Hooft for an iUuminating discussion on this. 
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below, consistent with effective bosonic string model predictions. At short distances the 
data matches perturbation theory. Therefore c(f) crosses the asymptotic value at some 
intermediate point. For SU{2), this intermediate point seems to be around 2ro or 1 fermi. 
Here it is still not clear whether a string behaviour has set in. It is only at a still larger 
distance of about 2.75ro, that the data seems to be well described by the Arvis curve. 

An interesting question is what happens to these scales for SU{N) as N increases. It 
has been observed in [18| that the intermediate qq separation at which c(f) assumes its 



asymptotic value decreases with increasing at finite lattice spacing. For SU (3) the data 



on the coarsest lattice in ref. just about touches the Arvis curve at l.SrQ. However 
both these scales shift towards larger r as one approaches the continuum limit. 

Comparing the data with effective strings at higher orders may tell us if the scale of 
string formation is the one suggested by the Arvis curve or happens earlier. From 0(r~^) 
and beyond, effective string theories motivate parametrising the leading deviations to c(f) 
from the Arvis behaviour as 

Ac{r) = Ar-^)\BrJlY. (4.3) 



r 

From our data in the region 1.6 — 2.4 ro we obtain, as best fit values, A = 0.209(9), = 
—0.235(24). The corresponding terms of the Arvis potential are ^Arvis = —0.00725, i?Arvis = 
—0.00145. It is of interest to know what the predictions of effective string theories are for 
these. 

At intermediate distances over a wide region of r varying from 0.5ro to 2.8ro, the data 
is very well described by a formula of the type 

-^-(i-^ + i) (") 



with a = 0.444(4), b = -0.258(2) and n = 0.357(15). Existing SU{3) data in 3-d ]14| 
also admit a similar description with nearly same values of a and b, but with a different n. 
The curves are shown in Fig. ^. At the moment it is not at all clear if there is any theo- 
retical basis for such a description. However they certainly provide accurate interpolation 
formulae. 

5. Conclusions 

In this article we have looked at the continuum limit behaviour of the SU(2) flux tube at 
intermediate distances by measuring the static qq potential. Starting from a distance of 
about 0.1 fermi where the potential starts breaking away from 1-loop perturbation theory, 
we go to distances of about 1.4 fermi where the data begins to be well described by the 
Arvis potential. 

We look at the continuum limit of the string tension and find complete agreement 



between open and closed strings |27]. Our data on c(f) seems to approach the Liischer 
term from below as expected in bosonic string models. At distances below 0.15 fermi the 
data joins onto the perturbative values. A modified Lennard-Jones type empirical formula 
describes the data well in the intermediate region. 
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Figure 3: Curves of the type given in eqn. 4.4 describing the c(f) data in the intermediate region. 
The red + is the continuum limit 3 — d SU (2) data. The green x , the blue * and the magenta □ 
are 3 — dSU{3) data from with tq = 3.30, 4.83 and 6.71 respectively. The curves LJl, LJ2 and 
LJ3 are given by (a, 5, n) of (0.444, -0.258, 0.357) , (0.458, -0.289, 0.691) and (0.442, -0.287, 0.498) 
respectively. The perturbative curves marked 1-p, 2-p^ and 2-p^B are given by the first, the 



first two and all the three terms of the expression: 



graCp r 



with 



Ait ro 2 

A = 0.013162 and B — 0.001089. The green filled squares denote the smallest possible values of 
r/ro for which c(f) can be calculated on that particular lattice. 



Further directions of study include increasing the qq separation to confirm that the 
data indeed stays on the Arvis curve. For 5f7(3) and SU{5), it would be really interesting 
to push to the continuum limit to see if the behaviour seen in SU (2) holds in those cases 
and find the distance where the data meets the Arvis curve. 
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4777 


0.021596 


(5) 


11.96 


-0 


12409 


(51) 


13 


12 


4796 


0.021450 


(6) 


12.96 


-0 


12580 


(67) 


14 


13 


4812 


0.021333 


(4) 


13.96 


-0 


12696 


(60) 


15 


14 


4825 


0.021240 


(4) 


14.97 


-0 


12819 


(86) 


16 


15 


4837 


0.021163 


(5) 


15.97 


-0 


12875 


(113) 


17 


16 


4847 


0.021100 


(5) 


16.97 


-0 


13032 


(154) 


18 


17 


4856 


0.021047 


(5) 


17.97 


-0 


13076 


(234) 


19 


18 


4864 


0.021002 


(6) 











Table 5: Force and c(f) (/3 10) 



r 




r 


fir) 




r 


c(f) 




3 


2 


3790 


0.024521 


(3) 


2.808 


-0.05095 


(1) 


4 


3 


4071 


0.019917 


(4) 


3.838 


-0.06530 


(3) 


5 


4 


4322 


0.017607 


(5) 


4.876 


-0.07714 


(5) 


6 


5 


4481 


0.016276 


(5) 


5.903 


-0.08677 


(10) 


7 


6 


4579 


0.015432 


(6) 


6.920 


-0.09460 


(15) 


8 


7 


4645 


0.014861 


(7) 


7.932 


-0.10105 


(23) 


9 


8 


4692 


0.014453 


(4) 


8.941 


-0.10640 


(17) 


10 


9 


4727 


0.014155 


(5) 


9.948 


-0.11081 


(25) 


11 


10 


4755 


0.013930 


(5) 


10.95 


-0.11446 


(35) 


12 


11 


4777 


0.013756 


(5) 


11.96 


-0.11744 


(47) 


13 


12 


4796 


0.013618 


(6) 


12.96 


-0.12005 


(55) 


14 


13 


4812 


0.013508 


(6) 









Table 6: Force and c(f) (/? = 12.5) 
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